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Abstract 



We construct a new exactly solvable supersymmetric spin chain related to the BCn 
extended root system, which includes as a particular case the BCn version of the 
Polychronakos-Frahm spin chain. We also introduce a supersymmetric spin dynam- 
ical model of Calogero type which yields the new chain in the large coupling limit. 
A ■ This connection is exploited to derive two different closed-form expressions for the 

O I chain's partition function by means of Polychronakos's freezing trick. We establish 

a boson-fermion duality relation for the new chain's spectrum, which is in fact valid 
for a large class of (not necessarily integrable) spin chains of BCn type. The exact 
Kj^ , expressions for the partition function are also used to study the chain's spectrum 

Hd I as a whole, showing that the level density is normally distributed even for a moder- 

O^ ■ ately large number of particles. We also determine a simple analytic approximation 

to the distribution of normalized spacings between consecutive levels which fits the 



numerical data with remarkable accuracy. Our results provide further evidence that 



f^ ' spin chains of Haldane-Shastry type are exceptional integrable models, in the sense 

OO . that their spacings distribution is not Poissonian, as posited by the Berry-Tabor 

^■"■^ ' conjecture for "generic" quantum integrable systems. 
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1 Introduction 



In the last few years exactly solvable (or integrable) supersymmetric spin 
chains and their associated dynamical models have been the subject of ex- 
tensive research in connection with different topics of current interest, such 
as the theory of strongly correlated systems p!p] or the AdS-CFT correspon- 
dence [3] . Among these models, the supersymmetric versions of the celebrated 
Haldane-Shastry chain [il[5] and its rational counterpart proposed by Poly- 
chronakos [6] and Frahm [7] occupy a distinguished position, due to the rich 
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mathematical structures at the heart of their highly solvable character. The 
Haldane-Shastry (HS) chain was introduced in an attempt to construct a 
simple one-dimensional model whose ground state coincided with Gutzwiller's 
variational wave function for the Hubbard model in the limit of large on-site 
interaction [8|9llT0] . It can also be obtained in this limit from the Hubbard 
model with long-range hopping studied in Ref. [11], in the half-filling regime. 
The original HS chain describes A^ spin 1/2 particles equally spaced on a cir- 
cle, with pairwise interactions inversely proportional to the square of the chord 
distance. Its Hamiltonian can be written as 

where J, = ^{af,af,a^), af is a Pauli matrix at site i, and the summa- 
tion indices run from 1 to A^ (as always hereafter, unless otherwise stated). 
A natural generalization of this chain to su{m) spin P^PH] is obtained by 
taking Jj = (J/, . . . , J™ ~^), where { Jf } are the generators of the fundamen- 
tal representation of su(m) at site i. In fact, with the usual normalization 
tr(J^Jj^) = i^'^T we have 

^ 2 V ^ m 

where Sij is the operator permuting the i-th and j-th spins, so that the su(m) 
spin Hamiltonian ([1]) is a linear combination of spin permutation operators. 

The spectrum of the spin 1/2 HS Hamiltonian was numerically analyzed in the 
original papers of Haldane and Shastry. The more general su(m) spin case was 
studied in a subsequent publication by Haldane et al. [13], who empirically 
found a complete description of the spectrum and explained its high degener- 
acy by an underlying y{slm) Yangian symmetry. These results were rigorously 
established in Ref. p3J by constructing a transfer matrix using the Dunkl 
operators p^fTT] of the (trigonometric) Sutherland dynamical model [T8|IT9] . 
Although this approach yields an explicit formula for the energies in terms 
of the so-called motifs, the computation of their corresponding degeneracies 
becomes quite cumbersome when m > 2. This is probably the reason why the 
partition function of this chain was computed only very recently pO], using a 
procedure known as Polychronakos's freezing trick which does not rely on the 
explicit knowledge of the spectrum [6f2T] . 

The key idea behind the freezing trick is exploiting the connection between 
the HS chain and the Sutherland spin dynamical model in the strong coupling 
limit. Indeed, in this limit the particles in the latter model "freeze" at the co- 
ordinates of the (unique) equilibrium of the scalar part of the potential, which 
are exactly the chain sites. Thus, in this limit the dynamical and the spin de- 
grees of freedom decouple, so that the spectrum of the spin Sutherland model 



becomes approximately the sum of the spectra of the scalar Sutherland model 
and the HS chain. This observation yields an explicit formula for the partition 
function of the HS spin chain as the strong coupling limit of the quotient of 
the partition functions of the spin and the scalar Sutherland models, both of 
which can be computed in this limit. In fact, this method was first applied 
by Polychronakos to evaluate the partition function of the su(m) spin chain 
associated with the spin Calogero (rational) model [22|23] . usually known in 
the literature as the Polychronakos-Frahm (PF) chain, whose Hamiltonian (in 
the antiferromagnetic case) is given by 

n = J:ic^-Q-^^+s^,)■ (2) 

i<j 

The sites Q of this chain are the coordinates of the equilibrium of the scalar 
part of the Calogero potential, which coincide with the zeros of the Hermite 
polynomial of degree A^ 



Both the HS and the PF chains admit a natural su{m\n) supersymmetric 
version, in which there are m bosonic and n fermionic degrees of freedom at 
each site [2^26|27|28|29j . In practice, the Hamiltonians of these chains are 
obtained by replacing the spin permutation operators Sij by suitable super- 
symmetric generalizations thereof (cf. Eq. (jl]) below). The partition functions 
of the su{m\n) PF and HS chains have also been computed in closed form us- 
ing the freezing trick method outlined above [25)127] . These partition functions 
were then used to establish a remarkable boson-fermion duality relation for 
the spectrum of each of these chains 



The original HS and PF chains (and their supersymmetric versions) are as- 
sociated with the Am-1 root system, since the interaction between any two 
spins depends only on the difference of their site coordinates. It is also of in- 
terest to consider generalizations of these chains to the BCi\f extended root 
system, in which the spins interact not only among themselves but also with 
their mirror images with respect to a reflecting end located at the origin. In 
the non-supersymmetric case, the HS chain of BCn type was first discussed 
in Ref. [30], and its partition function was evaluated by means of the freez- 
ing trick in Ref. [SI]. Likewise, the partition function of the su(m) PF chain 
of -BCat type introduced in [52] was recently computed in closed form by the 
authors [^. To the best of our knowledge, however, the supersymmetric coun- 
terparts of both the HS and PF chains of BCn type have not been studied 
so far. One of the main purposes of this paper is precisely to fill this gap in 
the case of the PF chain. More precisely, we shall apply the freezing trick to 
compute the partition function of this chain in closed form, relating it to the 
corresponding partition function of the su{m\n) PF chain of Ajv_i type. We 
shall also extend the boson-fermion duality mentioned above to a large class 
of supersymmetric spin chains of BCn type which includes the PF chain as a 
particular case. 



The spectra of the non-supersymnietric spin chains of Haldane-Shastry type 
whose partition function has been exphcitly computed share some statisti- 
cal properties that we shall now discuss. In the first place, the density of 
energy levels follows the Gaussian law with remarkable accuracy even for a 
moderately large number of particles [20|31|33IIM] . Secondly (except for the 
HS chain of BCn type, which has not yet been studied in this respect), the 
density of spacingsp(s) between consecutive (normalized) levels follows a char- 
acteristic distribution essentially different from Poisson's law p{s) = e"**. This 
fact is somewhat surprising, in view of the integrability of these chains 



and the long-standing conjecture of Berry and Tabor [35J, according to which 
the spacings distribution of a "generic" quantum integrable system should be 
Poissonian. The validity of this conjecture has been verified for a number of 
important integrable systems, including the Heisenberg chain, the t-J model, 
the Hubbard model [3S], and the chiral Potts model [37j. On the other hand, 
the non- Poissonian character of the spacings distribution has recently been ob- 
served for the su(m|ra) HS chain of A^^i type [27], whose density of spacings 
is qualitatively similar to that of the non-supersymmetric HS and PF chains. 
Another important aim of this paper is the study of the level density and the 
spacings distribution for the su(m|?T,) PF chain of BC^ type. Our main con- 
clusion in this respect is that this chain behaves in essentially the same way 
as other chains of Haldane-Shastry type previously studied. In other words, 
when the number of spins is sufficiently large the level density follows with 
great accuracy the Gaussian distribution, and the spacings distribution also 
deviates substantially from Poisson's law. As a matter of fact, we shall see 
that the spacings density is well approximated by the same "square root of a 
logarithm" distribution as the original HS and PF chains, and the su(m) PF 
chain of BCn type [33|IM] . This lends further support to the fact that spin 
chains of Haldane-Shastry type are exceptional integrable systems from the 
point of view of the Berry-Tabor conjecture. 

The paper is organized as follows. In Section [2] we recall the definition of 
the super symmetric spin permutation operators [25], and introduce new su- 
persymmetric spin reversal operators. With the help of these operators, we 
construct the su{m\n) PF chain of BC^ type and explain how the freezing 
trick can be used to compute its partition function from that of its associ- 
ated spin dynamical model. In Section [3] we compute the spectrum of the 
latter model, which is then used in Section H] to evaluate the chain's partition 
function. Using the grand canonical partition function of the spin dynamical 
model, in Section [5] we obtain an alternative expression for the chain's par- 
tition function in terms of that of the su(m|n) PF chain of A^-i type. This 
expression turns out to be particularly efficient for the computation of the 
chain's spectrum for a relatively large number of spins. Section [6] is devoted to 
generalizing to spin chains of BCn type the boson-fermion duality uncovered 
in Refs. [25|28] . In the last section we use the explicit formulas for the partition 
function obtained in Section \5\ to study the asymptotic behavior of the level 



density and the spacings distribution when the number of particles is very 
large. The paper also includes two appendices, the first one being a collection 
of several useful g-number identities, and the second a detailed derivation of 
the exact formulas for the mean and standard deviation of the chain's energy 
used in Section [71 



2 The model 



Let us begin by defining the Hilbert space of the internal degrees of freedom 
and the action therein of the su{m\n) permutation and spin reversal operators 
Slj and 5*"'. There are two species of particles, "bosons" and "fermions" 
(the reason for this terminology shall be apparent later), with m and n denot- 
ing respectively the number of bosonic and fermionic degrees of freedom. More 
formally, the internal Hilbert space S*^™'"^ ^ (^(^m+n^^N jg spanned by states 
of the form |si) ® ■ ■ ■ (8> \siy) = |si, . . . , sjy) = |s) , where Sj is a two-vector 
(slsj) with 

, „ , , 1 0, bosons 
i) s^ = 7r(sJ = < 

1, fermions 



11 s e 



{-^,-^ + 1,...,^}, 7r{s,) = l 



We shall say that s] is the value of the "spin" Sj, while s^ = n{si) is its type 
(bosonic or fermionic). As usual, the inner product on E'^'"'") is defined in such 
a way that the set of vectors |s) forms an orthonormal basis. 

The action of the spin reversal operators S"" on this basis is given by 

Sf\si, ...,Si,...,SN)= p{si) \si, ...,S~,...,Sn), (3) 

where s" = {-s},sf), 

Je, 7i{si) =0 

[e , TT{Si) = 1, 

and e, e' = ±1 are two independent signs. In other words, 5"" reverses the 
value of the i-th spin without affecting its type, and multiplies the state by 
the sign p(sj). 

The definition of the spin permutation operators S^j = S^^ is a bit more 
involved, namely (assuming that i < j) 

Oj^j I . . . , Sj, . . . , Sj, . . . / (^ L) I . . . , Sj, . . . , Sj, . . . / , (^4 J 



where 



and 



CtijiS 



(S) = n{Si)7l{Sj) + (7r(Si) + 7T{Sj))fij{ 



k 



i-1 

k=i+l 



is the number of fermions in between the z-th and j-th spins. In other words, 
the sign (— 1)"«j(^) is —1 when either both Sj and Sj are fermionic, or Si and 
Sj are spins of different type with an odd number of fermionic spins between 
them. An equivalent way of defining the action of Sj^ is by requiring that 

Slfj^i introduce a factor of —1 precisely when Sj and Sj+i are fermionic. This 
clearly implies the previous (apparently more general) rule, since an arbitrary 
permutation is a product of permutations of consecutive elements. 

The Hamiltonian of the supersymmetric sn{m\n) PF chain of BCn type is 
defined bjo 



n 



(m|n) 



E 



Q{m\n) -. 

^ + 



s. 



(m|n) - 



Lte-e,) 



te + Q' 



+ /?E' 



sr 



(5) 



where S, 



(m|n) 



■,(m|n) 



.jj — 5" 5" S'jj and /3 > 0. The chain sites ^j are the same as 
those of the ordinary PF chain of BCn type studied in Ref . [33] , namely the 
coordinates of the (unique) equilibrium of the scalar potential 



f/(x) = E 



1 



+ 



1 



yXi •^j) \^i \ Xj) 



+ E 



xf 4 



(6) 



{x = (xi, . . . , xjy) I < Xi < ■ ■ ■ < xat}. It can be 
/2yi, where i/i is the i-th zero of the Laguerre polynomial 



in the open set C 

shown [32] that C,i = 

L^ . Note also that the convention for the sign of the SUSY spin permutation 

operators in the definition of H}^ is the opposite to that of Ref. [27], but 

(as we shall see in a moment) our choice is more consistent with the use of 

the names "bosonic" and "fermionic" for the two species of spins. Since 



S, 



(m|0) 



1-3 



^J1 5 



Q(0|n) 



-s, 



ij ) 



sr = eS,, 



(7) 



where Si is the operator reversing the i-th spin, the purely bosonic (resp. purely 
fermionic) spin chain 7i(''"l*^) (resp. Ti^^'"^) coincides with the ordinary ferro- 
magnetic (resp. antiferromagnetic) PF chain of BCn type [33]. 



^ Note that the operators 5"^}™'"-' actually depend on e and e', although for simplicity 
we have chosen not to make this dependence explicit in the notation. 



The spin chain ([5]) is naturally related to the su(?Ti|ri) spin dynamical model 



jj{m\n) 



where xf. 
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a — S, 



(m|n) 
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s. 



(m\n) 



V 



X,- 



[X, 



I] J 



+ &E 



h-S^ 






4 



j3a and a > 1/2, and to its scalar counterpart 



Hf,, 



-A 






[X. 



-\2 



+ 



[X- 



K&-i)E4 + T 



i n 



Indeed, we have 



rr,{m.\n) 



TT{m\n) 



Hsc + a 7i„, 



(m|n) 



(9) 



(10) 



where Til^, is obtained from the spin chain Hamiltonian ([5]) by replacing the 

sites ^i by the coordinates Xj. Since 

i7H") = -A + a2f/ + 0(a), 



as a grows to infinity the particles tend to concentrate at the chain sites ^j, 
so that the dynamical and internal degrees of freedom decouple. By Eq. (fTOj) . 
in this limit the energies of the dynamical spin model ([8]) are approximately 
given by 



-^u — 



El" + a£j, a > 1 



(11) 

where E^" and Sj are any two energies of the scalar Hamiltonian ([H]) and 
the spin chain ([5]). The latter formula cannot be used directly to deduce the 
spectrum of TiJ^ from those of Hj^ and ifgc, since it is not known a priori 
which eigenvalues of H^^ and H^c combine to yield an eigenvalue of H^ . 
However, Eq. (fTTj) immediately leads to the exact "freezing trick" formula [2T] 



^{m\n) 



(T) 



lim 

a— >oo 



zi7i"Vr) 



'12) 



{rn\n) y{m\n) 



Am\n) 



relating the partition functions Z^^J''"'', Z'^l^^''''' and Zgc of H^J^^'' , H^'^'"'' and 
Hsc, respectively. The partition function Z^^ of the scalar model (^ was com- 

Am\n) 



puted in Ref. [33]. We shall see in the next section that the spectrum of H^^, 
is quite simple, a fact that shall be exploited in Section H] to explicitly com- 
pute its partition function. Using Eq. (fT2|) we shall then obtain the partition 
function Z^'^' in closed form. 



3 Spectrum of the spin dynamical model 



The first step in the evaluation of the partition function of the spin chain ([5]) by 
means of the freezing trick relation (iT2l) consists in determining the spectrum 



of the spin dynamical model ([8]). To this end, we introduce the auxiliary 
"exchange Hamiltonian" 



h'=-a+j: 






i^j l\-^i3J y-^'i-j) 



-T.^Ab-K.) + jr\ (13) 



where Kij is the operator permuting the z-th and j-th spatial coordinates, Ki 
reverses the sign of the z-th coordinate, and Kij = KiKjKij. The rationale for 
introducing this operator is the fact that its restriction to states symmetric 
under the action of the operators 

nS;!"^ = K^.slf'^^ , W' =K.ST\ z, J = 1, . . . , iV , (14) 

coincides with that of the Hamiltonian H^"^ , as we shall next discuss. To 
begin with, note that all three sets of operators {Ki,Kij}, {S^'"' , S^j } and 

{n"',n^J* } generate a realization of the Weyl group of BCn type, namely 
they satisfy the commutation relations 



Kij = 1, KijKjk = KikKij = KjkKik, K^Km = K^iKij, 

-2 



(15) 
Kt = 1, KiK^ = K,Ki , Ki.Kj, = K^K,,, K,,K, = K,K,, , 



(where i,j,k,l are distinct) and similarly for the other two sets. We shall 
denote by Aj^ the projector on states totally symmetric under the action 
of the operators {H"', I[\Y }■ For instance, if A^ = 3 we have 

A;-i") = l(i + nr')(i + nr')(i + nr) 



1 /. TT{m\n) T-r(m|n) TT(m|n) ■pr{m|n)-pr(m|n) ■pr{m|n)-pr(m|n) 



3! 

In particular, A^™l°^ (resp. A^°l"^) projects onto states totally symmetric (resp. 
antisymmetric) under particle permutations and with parity e with respect to 
simultaneous sign reversals of spatial coordinates and spins. From the defini- 
tion of the total symmetrizer A^™ it follows that 

yr{m\n) .{m\n) _ . (m|n)-p.(m|n) _ . (m|n) Tree' a (ni|n) _ \{m\n)yree' _ a (m|n) 

SO that 

^(™|n)^(rn|„) ^ ^^^^Mn) ^ sf A^J^^'''^ = K,A^^^'^ . (16) 

Hence 

which, as we shall see, is precisely the relation needed to deduce the spectrum 
of H^J^^""^ from that of H' . 



Indeed, recall [33] that the operator H' is represented by an upper triangular 
matrix in the (non-orthonormal) basis with elements 



0k = e-t^^ni^^MIK-^?| -114% (18) 

i i<j i 

where 

k={h,...,kN)EiNU{0}f = X. 
In other words, 

if '0k = -E^^k + Yl Ck'k0k' , 

|k'|<|k| 
where the eigenvalues are given by 

^k = a\k\ + Eo, 

with |k| = J2i h and Eq = Na{b + a{N -l) + fj. 

Due to the impenetrable nature of the singularities of the Hamiltonian ([H]), 
its Hilbert space is the set Ll(C) E^"*'"-* of spin functions square integrable 
on the open set C and vanishing sufficiently fast on the singular hyperplanes 
Xi±Xj = 0. However, it can be shown that H^l^ is equivalent to its symmetric 
extension (under nj™'"^ and Uf) to the space Ll(R^) ® s(™l"). The states 

^k,s = A;?"X0k(x)|s)), (19) 

form a basis of this Hilbert space provided that the following conditions hold: 

(i) ki ^ ■ ■ ■ ^ kN 
(ii) If ki = kj and i < j, then: 
(ii.a) 7r(si) ^ tt{sj) 

(ii.b) sj ^ sj + 7r(sj) whenever n{si) = tt{sj) 
(iii) si ^ for all i, and sj > if {-lf^p{si) = -1 

The last two conditions, which differ from the analogous ones for the su(?7i) 
case [33], deserve further remark. Condition (ii.a) simply states that acting 
with appropriate generalized spin permutation operators S^j we can first 
reorder the spins so that within each sector of s (spin components correspond- 
ing to a maximal sequence of equal components of k) bosons precede fermions. 
Furthermore (condition (ii.b)), the bosonic spin components can be arranged 
in non-increasing order, while for the fermionic ones we can enforce strictly de- 
creasing order since in this case the state vanishes if Sj = Sj by antisymmetry. 
(This is, by the way, the justification for the names "bosonic" and "fermionic" 
used for the two species of spins.) Finally, condition (iii) stems from the fact 
that the spin functions flT9l) must be even under Hf , while if s} = we have 

nr>k,s = Ai?"^(nr'0k(x)|s)) = (-i)^v(^.)V^k,s, 



so that (— l)^*p(si) = 1. A direct consequence of this condition is that the 
number of bosonic spin values of Sj compatible with rule (iii) is given by 

y , m even 

m{kj) = {^, m odd, k. even (20) 



m—€ 



m odd, fci odd . 



2 ' yj^yj-, 'i^j 

Similarly, the number n{kj) of fermionic spin values of Sj compatible with the 
third rule is given by the previous expression, with m replaced by n and e 
bye'. 

The Hamiltonian H^ is represented by an upper triangular matrix in the 
basis (IT^ . ordered according to the degree |k|. Indeed, taking into account 
that H' commutes with A),™ , by Eq. fITTI) we have 

|k'|<|k| 

Hence the eigenvalues of H^^^ are given by 

Ek,s = a|k| + Eo , (21) 

where k and s satisfy conditions (i)-(iii) above. Note that, although the value 
of i?k,s is independent of s, the degeneracy of each level clearly depends on 
the spin through the latter conditions. 

From Eqs. ([Tjj-Q it follows that the Hamiltonian H^c of the scalar Calogero 
model is simply H^ , so that its spectrum is also given by the RHS of 
Eq. (1^ . with all ki even by condition (iii). We can thus subtract the ground 
energy Eq from the spectra of both H^c and H^^ without altering the parti- 
tion function of the chain 1^. With this proviso, the energies of both models 
are proportional to the coupling constant a, which implies that Zsc(aT) and 
Z^^ (aT) are independent of a. Hence, in the calculation of the partition 
function of the chain ([5]) we can take without loss of generality Eq = and 
a = 1, and write the freezing trick formula flT2|) simply as 

4™'"^(T) = ^|-^. (22) 

4 Evaluation of the partition function 



We are now ready to compute the partition function of the spin chain ([5]) 
using the freezing trick formula (1221) . As for the ordinary (purely bosonic or 
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fermionic) chain, the evaluation of the partition function Zj^ of the spin 
dynamical model ([8]) depends on the parity of the integers m and ra, so that 
one has to consider four different cases. In the calculations that follow, we 
shall use the fact that the number of ways into which we can arrange the spin 
components Sj in a sector of s corresponding to a certain component «; of k 
repeated v times is given by 

Clearly, this number depends on k only through its parity, cf. fl2Up and its 
fermionic analogue. Note also that when m{K) = the second binomial coef- 
ficient in fl23|) should be interpreted as 6i^. 



Case 1: m,n even 

Recall, to begin with, that after setting Eq = and a = 1 the partition 
function of the scalar model ([9]) reads [33] 

ZsciT)= y: ?'"' = na - 9'r \ g=e-v(^BT)_ (24) 

By Eq. (|2T1) . the partition function of the Hamiltonian ([8]) can be written as 

4™'"H^)= E ^k^l, (25) 

where the spin degeneracy factor t/k is the number of spin states |s) satisfying 
conditions (ii) and (iii) for a given multi-index k. Since both m and n are even, 
by Eq. (!20l) and its fermionic analogue we have rn{kj) = m/2 and n{kj) = n/2. 
If 



k = (Ki, . . . ,Ki, . . . ,Kr, . . . , Kj.), Ki > ■ ■ ■ > Kr ^ 0, (26) 

from Eq. (!23|) it follows that 



d^=lld{u^)=d{iy), (27) 



where v = {ui, . . . ,h'r) and 

Since J2l=i ^i = -^5 the multi-index i/ is an element of the set Vn of partitions 
of A'" with order taken into account. Inserting Eqs. (I26l) - (!28l) into (l25l) we 
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obtain 

r 

The RHS of the previous equation can be evaluated as in [331 Eq. (24)], with 
the result 

ZiMn) ^ ^-N ^ ^(^) Jl ^^ ^ (29) 

where £(i>') = r is the number of components of the multi-index u and 

E 



iv, = E^.- 



From Eqs. (l22l) . (12^ and (!29|) it follows that in this case the partition function 
of the su(m|ra) PF chain of BCn type is given by 

t(u)-l 

J2 Nj N-i{u) 

where the positive integers N'j are defined by 

{iv(,...,iv;_,(,)} = {i,...,iv-i}-{iVi,...,iv,(.)_i}. 

Note that, as for the purely bosonic/fermionic model, in this case the partition 
function does not depend on e, e'. We shall therefore drop from now on the 
subscripts e, e' from Z^^ when both m and n are even. The first product in 
Eq. (l30l) is simply the partition function of the su(2|0) chain ([5]) (cf. Eq. (34) 
in Ref. [33]), whereas the remaining sum coincides with the partition function 

Zj^Y^^ of the su(f If) PF chain of type A^-i computed in [29]. We thus 
obtain the remarkable relation 

^H") = 2(210)^^12)^ m,ne2Mo. (31) 



Recall [25] that the partition function of the su(k\l) PF chain of type A is 
explicitly given by 

2"">- r /!."■'"'-" M" (32) 

where the notation {q)k is defined in Eq. ( JA.ll) of Appendix Rl For instance, 
from Eq. ([32]) with k = l = l and Eqs. flATB - fTOj) one obtains 



N 

^(A) - 2^ 

M=0 



AT' 
M 



N-l 



gUJ= n(i + ?' 

g i=0 
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so that 

N-l 

^(2|2) ^ ^(210) ^(i|i) ^ 2(1 + q^) n (1 + gT . (33) 



Case ^; m oo?c?, n even 



In this case, by Eq. (120|) and its fermionic counterpart we have n(fcj) = n/2, 
while m{kj) depends on the parity of kj. Hence it is convenient to modify 
condition (i) by first grouping separately the components of k with the same 
parity and then ordering separately the even and odd components. In other 
words, we shall now write k = (ke,ko), where 



kg — I 2/?!, . . . , ^Kl, . . . , Zl^g, . . . , Zl^s j , 



ko = ( 2Ks+i + 1, . . . , 2Ks+i + l,...,2K.r + l,...,2K,r + l 

and 

Ki > ■ ■ • > Ks ^ 0, Ks+i > ■ ■ ■ > Kr ^ . 

ByEq. l^,m{2Kj) = (m + e)/2forj = 1, ..., s, while m(2fi:j + l) = (m-e)/2 



for j = s + 1, . . . ,r. From Eq. fl23l) it follows that the spin degeneracy factor 
dk is now given by 

j=l j=s+l 

where 

Equation ( !25l) for the partition function of the dynamical spin model ([H]) then 
becomes 

s r 

I^ePjV S = Kl>--->Ks>0 

where the RHS can be evaluated as in [33, Eq. (29)]. We thus obtain 

^H.)^ E E^:(-)^-^"""^^nY3^- n 7^^, (36) 

where 

j 



i=s+l 
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Substituting Eqs. (^^ and (136|) into (l22l) we arrive at the following explicit 
expression for the partition function of the spin chain ([5]) for odd m and even n: 

i{i^) s 2N, t{^) 27V, 

Note that in this case the partition function depends on e but not on e'. 

As is the case for any finite system with integer energies, the partition function 
Z^^ should be a polynomial in q. We shall now prove this fact by simplifying 
Eq. (137|) with the help of the identities in Appendix |Al To this end, we define 
the two sets of integers 

[n[, ..., N'^^_,} = {1, . . . , AT, - 1} - {aTi, . . . , Ar,_i} , 
{iV^,_.+i,...,iV^_,H} = {l,...,Ar-Ar,-l}-{iV,+i,...,iV,(.)_i}, 

in terms of which we can rewrite Eq. (1371) as 



f(,.)-i 



, , . ^H Af-JVs+2 V Afj+2 V Nj 






Ns-s N-e{u) 

i=l j=Ns-s+l 



From the discussion in Appendix |A] (cf. Eqs. (lA.6p - (]A.7P ). it follows that the 

is an even polynomial of degree 2Ns{N — Ng) in q. 



binomial coefficient 



Thus the RHS of Eq. (155]) is also a polynomial in q, as expected. 
Case 3: m even, n odd 

The evaluation of the partition function z^^ is performed as in the previous 
case, the only difference being that the spin degeneracy factor dl{u) in Eq. (1551) 
should be replaced by 

<,H = n4'(^.)- n d-Ai^j), (39) 

J = l j = s + l 

where now 

In particular, the partition function Z^^ does not depend on e in this case. 



14 



Case 4: m^n odd 

In this case both rn{kj) and n{kj) depend on the parity of kj, so that the 
computation of the partition function proceeds as in Case 2. The final result 
is still Eq. ( 155]) . with the spin degeneracy factor dl{u) replaced by 

dl,H = f[d,A^,)- n d^^M^,), (41) 

j=l j=s+l 

where 

Thus in this case the partition function depends on both e and e'. 



5 Relation with the spin chain of Ajv_i type 



We have seen in the previous section that when the integers m and n are 
both even, the partition function Z^"^'^^ of the chain ([5]) is related in a simple 
way to that of the su(y||) PF chain of A^^i type, cf. Eq. (I3T1) . The purpose 
of this section is to derive an analogous relation for odd values of m and/or 
n. To this end, we shall deduce an alternative expression for the partition 
function by means of the grand canonical partition function of the spin dy- 
namical model (IH]), following an approach similar to that of Refs. [2Tp5] for 
the ordinary and supersymmetric PF chains of A^^i type. 

In order to compute the grand canonical partition function, it is convenient 
to replace conditions (i)-(iii) in Section [3] ordering the labels k, s of the basis 
states (TT9|) by the following equivalent set of rules: 

(a) li i < j then 7r(sj) ^ vi"(sj) 

(b) li i < j and 7r(sj) = 7r(sj), then sj ^ s] 

(c) If z < j and Si = Sj, then ki ^ kj + 7r(sj) 

(d) Si ^ 0, and (-1)'=' = p{si) if sj = 0. 

By Eq. ( 1211) (setting, as before, Eq = and a = 1) the grand canonical 
partition function Zj^ of the spin dynamical model (jHD is given by 

oo 

z;?"^(r,/i) = EE5'""''"> (43) 

N=0 k,s 

where n is the chemical potential and the inner sum runs over all values of 
k G Nq and s compatible with conditions (a)-(d). As in the previous section. 
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the computation of Z^^ depends on the parity of m and n. Since the case 
of even m and n has aheady been dealt with above, we shall start with the 
case of rra odd and n even: 



Case 2: m odd, n even 



Calling m = {m — l)/2 and n = n/2, by conditions (a), (b), and (d) we can 
write 



l^fh + l 



^TTi + n 



0,... 



m—1 



m-1 1 



n-1 



n-1 



'•••'"'•••' 2 '•••' 2 52'''''2''''' 2 '■■■' 2 /' 



(44) 



where z/j G Nq and I]S)" ^i — ^- With a slight abuse of notation, we have 
suppressed the second component s^ of each spin Sj, since clearly s^ = for 
the first i^o + - ■ ■ + i^m spins and s^ = 1 for the remaining ones. By condition (c), 
the multi-index k is of the form 



k= 2/.^ +p,...,2<+p, /.!,...,<, ...,/.r",---, C7. 



where 



P 



> 



(45) 



(46) 



and the nonnegative integers k) satisfy 



t^\^ 



K\ > 



^< 



> K 



i = 0, . . . ,m 

i = m + 1, . . . ,m + "h . 



By Eq. (1431) . the grand canonical partition function of the chain ([5]) is given 
by 



"0 



Z{m\n) 



E E E ^'^^ n 5^=' 

l/Qv.l'A+n^O K»J;...J;4, 4>...>4. '=1 



^— 0,...,m j='/Ti+l,...,m+n 



E (!'=' 



n E <i 



E Kl-Ui^l 



1/0^0 



i^,^0 



X 



rh+h 



n E ^ 



j^ 

r = l 



Y^ ni-Ujfi 



The first factor in the last equality is clearly the grand canonical partition 
function Z^^'^-' of the purely bosonic spinless dynamical model ([8]) (which co- 
incides with the scalar model (19]) for e = 1). On the other hand, in Ref. [25] it 



16 



was shown that 

1^ = Kl>---^Ky !/=0 Kl>-->Ki/ 

where Z|^-, (resp. Z| A ) denotes the grand canonical partition function of the 
purely bosonic (resp. fermionic) spinless Calogero model of type A. We thus 
have 

Using again the results in Ref. [25j, we can write the previous formula as 

2Hn)_Z(i|o)2(yit)^ (47) 

where Z/^? ^ denotes the grand canonical partition function of the su(^^^ 1 1) 
spin Calogero model of type A. On the other hand, the grand canonical par- 
tition function Z{T,fi) can be expressed in terms of the A^-particle canonical 
partition function Zi\i{T) as 

oo 

Z(T,/i)=5:y^Z^(T), 

Af=0 

where y = q^^ and Zq{T) = 1. Equating the coefficient of y^ in both sides of 
Eq. (H7|) we readily obtain 

r^{m\n) _ ••r-^ 7(110)7(^^15) 
^ee',N — Z^ ^e,M ^{A),N-M ' 
M=0 

/ 771—1 I n \ 

where Z,^^j^^j^ is the partition function of the su(^^^ 1 1) spin Calogero model 
of Ajv_i type with N — M particles. By Eq. (!22|) . the partition function of the 
A^-particle su(m|?2) PF chain of BCn type is thus given by 

^ee',N - Z^ -^ ^{A),N-M ■ K^^J 

M=Q ^sc,N 

The partition function Z^ ^ is simply given by 

^+,M = Zsc,M = I ON ■. (49) 

where we have used (1241) . On the other hand, from conditions (i) and (iii) in 
Section [3] we have 

M 
, , ^ V(2fci + 1) „M 

^-M - Z^ <? - ? ^sc,M - / 2^1 • ^^^1 

ki>--^kM^O W )M 
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The partition function of the su(^^^^||) spin chain of type A with N — M 
particles is related to the partition function of its corresponding spin dynamical 
model by 

/ t-n — ^ 1 I '»-> \ / 'm — 1 I 'Ki s 

(51) 



/ m— 1 \ n \ / m — 1 \ n \ 

C^y 2 I 2 ' I „\ ryK 2 \2> 

^{A),N-M — y(l)N-M^^A),N~M ■ 

Inserting Eqs. (I49l) -( I5T1) into ( l48l) we finally obtain 



^{m\n) 



N 



(q' 



ee',N — Z^ y 
M=0 

N N 

i(l-€)Af 



z. 



f m — 1 I ri_\ 



Af=0 



(g2)M(g)iv-Af ^(A)'^-^'^ 

iV 



n (i+? 

j=M+l 



M 



57"- 2 I 2 •' 
^{A),N-M ' 



(52) 



where we have used the identity flA.2p . Alternatively, from the first equality 
in ( 152|) and Eq. (13T1) we can also derive the relation 



o-(m|n) 



AT 



E 1' 

M=0 



l(l-e)M 



AT' 
M 



Z 



(m— l|n) 

'N-M 



Case 3: m even, n odd 



This case is very similar to the previous one, the main difference being that 
Eq. (H71) now becomes 



y(m|n) _ y(0|l) y(Tl%^) 

which leads to the following expression for the chain's partition function: 



(53) 






(m\n) 



N 



N 7(011) 
M=0 ^sc,N 



(54) 



By conditions (i)-(iii) in Section [3l the partition function Z\_ \,l of the spinless 
fermionic model ([H]) with even parity is given by 



'(0|1) 

'+,M 



E (1 '=' 



2Y,{h+M-i) 

E ? '=' 



ki>--->kM^O 



h>-^lM>0 



Similarly, 



A/(M-l)v 



q 



M{M-1) 



sc,M 



ii': 



M 



■ (55) 



Ko|i) 



,Af^ 



ki>->kM>0 ^^ '^ 



(56) 



Inserting (155!) and (l56l) into Eq. (IMl) and proceeding as before, we obtain the 
following two expressions for the partition function of the chain ([5]) with even 
m and odd n: 



N 



N 



i=M+l 

'N 



M=0 

N 



N 
M 



/ m [ n — 1 \ 

5r\ 2 I 2 ' 
^(A),Ar-M 



J2 g*^(A'^-|(l+^')) 



M=0 



M 



z 



(m|n— 1) 



(57a) 
(57b) 



Jg- 



Case ^; m, n odd 



In this case both the bosonic and fermionic spin components can take the zero 
value, so that instead of Eq. P7|) we now have 



(l|0)y(0|l)^(^l^) 



Z(m|n) _ ry(l\V) ry\u\±j ry 
ee' — ^e ^e' ^{A) 

In particular, setting m = n = 1 in the previous this formula we have 

7(111) _ ymo) y{o\i) 

and therefore 

ry{m\n) _ ,7(111) ryl^^l^-j 

Proceeding again as in Case 2, we easily obtain 






N 7(111) ,__-,,„_-,. AT /2\ 



/j , ,, /m--l|n--l\ 

7 ^{A),N-M — Z^ /ox / N 



;r(l|l) ;r(VlV) 



(58) 

(59) 
(60) 

(61) 



As before, this expression can be simplified in two alternative ways, leading to 
the following remarkable identities for the partition function of the chain ([5]) 
with odd m and n: 



N N 

4^? = E n (1+9^ 

Af=0 i=M+l 



N 
M 



7(111) 7(^1^) 

^ee'.M ^{A),N-M 



N 

E 

Af=0 L 



M 



^(l|l) ^(m-lin-l) 



(62a) 
(62b) 



Let us finally evaluate the partition function Z_:^,' appearing in Eqs. (162!) . To 
begin with, from Eq. (l59l) for the grand canonical partition function and the 
freezing trick formula (!22|) we obtain 



N 7(110) 






|1) 



•7(111) _ Y^ ^e,N-M^e',M 



M=0 



^sc,7V 



TV 
y^ p{N-M)+M{M-l+p') 

M=0 



N 
M 



(63) 



Jg- 
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where p is given by (H6l) and p' = {1 — e')/2. We thus have 



N 



'{1|1) ^ V- „2( 

A/=0 



2ir=E9'^- 



M 



N-l 



2 11(1 + ?^ 



2i\ 



(64) 



jg- 



i=l 



where we have used Eq. (]A.3P with x = 1 and g replaced by g^. On the other 
hand, if e = e' = — 1 Eq. (IB^ becomes 



^a^) ^ ^iv ^ ^2( 



Af=0 



AT 
M 



„iV2(l|l) 



(65) 



Jg- 



Next, if e = — e' = 1 we have 



Z 



(i|i) 



N 



E^ 

Af=0 



Af2 



N 
M 



N-l 
g2 i=0 



(66) 



where we have used Eq. (IA.4p . Finally, when e = — e' = —1, Eqs. (|63|) and flA.SP 
yield 



N 

c 

M=0 



^ail) = gA^ ^ qM{M-2) 



N 
M 



N-l 

g^ n (1 + 9"'-') 



Jq- 



j=0 



Af-1 



g"'-'(l + 9)^11(1 + ?"'"')• (67) 



i=2 



6 Boson-fermion duality 

In Ref. [2B] it was uncovered a remarkable boson-fermion duality satisfied by 
the spectrum of any spin chain of Atv-i type with global su(m|ra) symmetry. 
We shall now extend these ideas to the PF chain of BCn type ([5]), although 
our results will in fact be valid for any chain of BCn type possessing global 
su(m|ra) symmetry. 

Let us start by defining the star operator * : S*^™"'"^ -^ ^(m|n) |^y 

|s)* = (-l)I^.^-(^>)|s), (68) 

and the exchange operator A"*^™!") : S(™l") -^ j^{n\m) |^-y 

A'Hn)|,^,...,,^) = |,;,...,4), ,^^(,i,i_,2^. (69) 



In other words, s'^ has the same value as Sj but opposite type. Following 
Ref. ^, we next define the operator ZY^™!") : E^™!") -^ £("1"^) as the com- 
position 

^(m|n) ^ ^{m\n) q * _ /^Q\ 
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Since, clearly, the operators * o * and A:'("l™) o ^("^1") are the identity in S(""l"), 
l/{Mn) is invertible and 

Indeed, if |s) G S*^"!™'^ we have 

Moreover, it is also easy to see that the star operator is self-adjoint and 

;f (m|n)t ^ ^(n\m)^ g^ ^j^g^^ ^(m|n)t ^ *o;f ("I™) = ^(^1")-^ ^ j.e., Z^C™!'^) is unitary. 

A straightforward computation [2B] shows that ui^l^) S^f""^ = -S\f"'^U^'^\''\ 
or equivalently 

The above considerations are enough for proving duality in the A]^ case. How- 
ever, in the BCj^ case we need also consider the behavior of the spin reversal 
operators S*" with respect to conjugation by W(™l"). For clarity's sake, in 
the following discussion we shall use the more precise notation s^ to 

denote the operator Sf : S^'"!") -^ S^"'!") . Noting that (s")' = [s[)~ and 
7r(sr) = 7r(si), if |s) G S(™l") we have 

and therefore 

Ti(m\n)~^c(ne'\me)i,[m\n) _ q(me\ne') _ _ n(.m,-e\n,-e') /r^^) 

With the help of Eqs. flTTl) and fl72l) one easily obtains 

where the last sum was evaluated in Ref. [33] • Since W^"*!") TY^r^ ^('"l") and 
^(n|mj ^^^ isospectral, from the previous equation we obtain the remarkable 
duality relation 

4'^H(g) = g^^zi7Ji(g-i). (73) 

When e' = e, the previous formula relates the spectra of two chains which 
differ in the exchange of bosons with fermions and the action of the spin 
reversal operators. On the other hand, when e' = — e Eq. (173!) is a genuine 
duality relation with respect to the exchange of bosons and fermions, i.e., it 



21 



relates the spectra of the chains 7i_^" and 7i_J7 • In addition, ii m = n 
Eq. (1751) establishes that the spectrum of the chain Tiill^' is invariant under 
the transformation £ t— > N'^ — S, i.e., it is symmetric about N'^/2. For instance, 
for m = 1 this property can be checked directly from Eqs. (!66l) - (!671) for the 
partition function. 

Similarly, if m and n are both even the partition function does not depend on 
e and e', so that Eq. ( ff3l) is again a genuine boson-fermion duality relation. 
For the same reason, ii m = n G 2No the spectrum of the chain TiJ^ is 
invariant under S i— > N'^ — S. As an example, if m = n = 2 this property can 
be readily verified using the exphcit formula Eq. ( l33l) . 

As mentioned above, a duality relation of the form (!73|) clearly holds for an 
arbitrary chain with global su{m\n) symmetry, with Hamiltonian 

where Cij, Cij and q are real constants. More precisely, from Eqs. (!7T!) and (172!) 
it now follows that 

so that Eq. ( 1731) becomes 

Z^V '(g) = g^Zi ' ^,(g ). 



7 Level density and spacings distribution 



Having derived several explicit formulas for the partition function of the 
spin chain (^, it is natural to inquire on the global properties of its spec- 
trum. In fact, the partition function (and hence the spectrum) can be com- 
puted in a very efficient way for relatively large values of N by expanding 
in powers of q Eqs. (IM]) , ([52]) , (IFTaj) . (l62iD and (I6l-([67D, together with 
Eq. (I32l) . with the help of a symbolic computation package. For instance, 
using Mathematica"'"'^ on a personal computer it takes about five seconds 

('21 1 '1 

to evaluate the partition function Z__ for A^ = 40 particles. 

In the first place, our calculations for a wide range of values of A^, m and n show 
that the energy levels of the chain ([5]) are equidistant, as already observed by 
us in the purely bosonic or fermionic case [33]. More precisely, we have found 
that the distance 6S between consecutive levels is equal to one, except for the 
su(l|l) chain with e = e', for which 6S = 2 hj Eqs. (1M|) and (1^ . 
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Secondly, our numerical calculations evidence that for sufficiently large A^ 
(A^ > 20) the normalized level density 



f{£) = {m + n)-''Y.di^i^~^i. 



(74) 



i=l 



where £i < ■ ■ ■ < £l are the distinct energy levels and di is the degeneracy of 
£i, can be approximated with great accuracy by the Gaussian law 



9{S) 



27ra 



e '^^^ 



(75) 



with parameters ^ and a given by the mean and standard deviation of the 
chain's spectrum. Since the energy levels are equidistant, this means that 



di 



6£{m + n 



\N 



9{£i 



(Ar>l) 



(76) 



As an illustration, in Fig. [T] we have plotted both sides of the latter equation in 
the case m = n = 1 and e' = — e = 1 for A^ = 15 and A^ = 30. The approximate 
Gaussian character of the level density is thus a property shared by all spin 
chains of HS type studied so far, both of type A and B [20|27f31f33f34] . 





Fig. 1. Plot of the Gaussian distribution (j75p (continuous red line) versus the LHS 
of Eq. (j76p (blue dots) in the case m = n = 1, e' = — e = 1, for A^ = 15 (left) 
and A^ = 30 (right). The root mean square errors (normalized to the mean) of the 
adjustments are 4.0 x 10~^ and 2.2 x 10~^, respectively. 



Since, by Eq. (!76|) . the level density is asymptotically given by the Gaussian 
law (!75|) . it is of interest to compute the mean energy and its variance as 
functions of A^, m, n, e and e', as has been done for other spin chains of 
HS type (including the purely fermionic version of the chain ^). Indeed, in 
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Appendix [B] we show that fi and a^ are given by 




(m + n) 



iV(iV-l)+('l- '^'" + '^" ')-, (77) 

V m + n J 2 



P^]^iAN^ + 6N-l) + ^^^NiN-l)iN-2) 

2 

{neprn - me'pn)N{N - 1) 



m + ny 

N ( fm — n^"^ 



4(m + nY \ \m + n 



(epm + e'Pn)' , (78) 



where Pm,Pn ^ {0, 1} are the parities of m and n, respectively. In particular, 
when N tends to infinity /i and a^ grow as A^^ and A^^, as for the ordinary 
(non-supersymnietric) PF chains of type A and B (cf. Refs. (SUES]). 

Another property of the chain's spectrum worth studying is the distribution of 
spacings between consecutive "unfolded" levels. In general [38], the unfolding 
of the levels £i of a spectrum is the mapping E^ \-^ rji = rj^Si), where r]{£) is 
the continuous part of the cumulative level density 

J-oo [m + n) ..^,^g 

It can be easily shown that the unfolded spectrum {r/j}^! is uniformly dis- 
tributed regardless of the initial level density, making it meaningful to compare 
different spectra. In our case, by the previous discussion we can take ri{S) as 
the cumulative Gaussian density fl75l) . namely 



J —CO Z 

One then defines the normalized spacings 



erf 



S — fi 
V2(J 



(79) 



Si = {r]i+i - r]i) / A , i = 1,...,L-1, 

where A = {rji — 771)/ (L — 1) is the mean spacing of the unfolded energies, so 
that {si}iSi has unit mean. 

As mentioned in the Introduction, the main motivation for studying the spac- 
ings density p{s) lies in the Berry-Tabor conjecture, which states that the 
distribution of spacings of a "generic" quantum integrable system should obey 
Poisson's law p{s) = e^'^. By contrast, for a chaotic system the spacings distri- 
bution is expected to follow Wigner's surmise p{s) = (vrs/2) exp(— 7rs^/4), as 
is approximately the case for the Gaussian ensembles in random matrix the- 
ory ^38j. A detailed study of the spacings distribution has been performed for 
many spin chains of HS type, namely the PF chains of types A and B [Mf33] . 
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the original (type A) HS chain [20JIM] . as well as the supersymmetric version of 
the latter chain [27]. The spacings distributions of all these chains, which turn 
out to be qualitatively very similar, differ essentially from both Wigner's and 
Poisson's distributions. More precisely, for the ordinary (non-supersymmetric) 
chains mentioned above we showed in our recent papers [33|3^ that the cu- 
mulative spacings distribution P{s) = Jo pW) ds' is approximately given by 




P(.)^l-^ llogR^h (80) 



where Smax is the maximum spacing. As a matter of fact, in Ref. [33] we proved 
that the previous approximation holds for any spectrum ^min = Si < ■ ■ ■ < 
£l = f^max, provided that the following conditions are satisfied: 

(i) The energies are equispaced, i.e., Si+i — Si = 6S ior i = 1, . . . , L — 1. 

(ii) The level density (normalized to unity) is approximately given by the 
Gaussian law fl75l). 



(iii) £^max -/",/"- ^min > CT- 

(iv) f^min and £^max are approximately symmetric with respect to fi, namely 

l^min T C-max ^f^l *^ ^max ^inin- 

Moreover, when these conditions are satisfied the maximum spacing can be 
estimated with great accuracy as 

S - £ ■ 

♦-Tnax '-'min /n-i \ 



271 a 

It should also be noted that Eq. ( IHOl) is valid only for spacings s ^ sq, where 
So is the unique zero of the RHS of this equation, namely 

So = Smax 6"?''^=''^ • (82) 

Since, by Eq. (IHT!) and condition (iii), we have (for instance) 

4 n^ax ^ g^2 ^ -^ ' 

it follows that sq <^ Smax- The validity of the above conditions for the chain ([5]) 
when N ^ 1 was checked in Ref. [33] in the purely fermionic case, which 
automatically implies their fulfillment in the purely bosonic case on account 
of the duality relation (175]) . For this reason, we shall assume in the rest of this 
section that m,n^ 1. We shall see below that in this case conditions (i)-(iii) 
are also satisfied when iV 3> 1, while condition (iv) only holds for m = n. 
It turns out, however, that when one drops condition (iv) the approximate 
formula (l80ll still holds, albeit in a slightly smaller range. 
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Indeed, it was proved in Ref. [33] that conditions (i)-(iii) imply that the spac- 
ings Si are approximately related to the energies £i by 

Si ^ Smax e ~^^^ . (83) 

Hence 

#(s. < ■•>) 1 



where 



#^mm ^ ^^ ^ ^- +#^+ ^ ^^ < ^1 



, (84) 



£± = /x±v^aWlog(^) (85) 

are the roots of the equation s = Smax e 2,72 ^ cf. Fig. [2l Using condition (i) 
to estimate the RHS of the latter equation we easily obtain 



Pis) 



1 



{L-l)5£ 
If 



max(£_ - £:min, 0) + max(£:max - ^+, 0) . (86) 



max ( Smax e 2.2 ^ ^^^^ e 2^2 



^7) 



then for s ^ s^i^ we have £_ ^ f^min and S^ ^ ^maxi see again Fig. [21 Hence, 
for s ^ Smin Eq. fl5^ becomes 



Pis) ~ ^_^ ^^ (^- - ^min + ^max " £+) ■ (88) 



Since (again by condition (i)) £^max — ^min = iL — l)5£, substituting Eq. (1851) 
into (l88l) and using (IHTl) we immediately arrive at the approximation (IHOl) 
for P(s). We have thus shown that conditions (i)-(iii) imply that Eq. flHOl) 
holds for s ^ Smin- Note, finally, that the minimum spacing Smin satisfies the 
inequalities 

Indeed, the second inequality is an immediate consequence of Eq. f lHTj) and 
condition (iii). As to the first one, by Eqs. flHTl) . fl82|) and flHTj) we have 



Q 2 V'-'iiiax '-'mmj 



log ( I = -— ^ (f^max - ^n 



< 



1 



2 min(/i - S^i^, ^max - /i) = log 



2 . /S 



8a^ 
where the equality holds if and only if S-^^n and ^max are symmetric about fj,. 
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Fig. 2. Plot of the RHS of Eq. (j83|) showing the points <Smirn "^maxi and the roots 8± 
corresponding to a spacing s ^ Smin- 

Let us next verify that the spectrum of the chain ([5]) satisfies conditions (i)- 
(iii) when A^ ^ 1 and m,n ^ 1. We have already seen at the beginning of 
this section that conditions (i) and (ii) hold for sufficiently large N . In order 
to check the validity of condition (iii), we need to compute in closed form 
the chain's maximum and minimum energies. By equation (ITT!) , the minimum 
energy is given by 

a-^oo a 

where -E^in ^^'^ -^min are the minimum energies of the scalar and spin dynami- 
cal models ([9]) and ([H]), respectively. From the discussion in Section [3] it follows 
that -Emin = -^0, so that (cf. Eq. ( 12T1) ) Smin is the minimum value of |k|, where 
k is any multi- index compatible with conditions (i)-(iii) in Section [3l If m > 1 
we can obviously take k = (and, for instance, s = f (s, 0), . . . , (s, 0) j, where 
s is any positive spin component), so that £^min = in this case. Similarly, 
when m = e = 1 the minimum energy is again zero, since the pair k = and 
s = f(0, 0), . . . , (0, 0) j satisfies the required conditions. However, if m = 1 and 
e = —1, condition (iii) implies that ki must be odd when s} = 0. Hence, the 
multi-index k yielding the minimum energy is of the form 




71(0) <N, 
rr(0) ^ A^, 



with Siain = max (A^ — n(0),Oj, where the notation n{kj) is defined by the 
fermionic analogue of Eq. (!20|) . (A compatible vector s is obtained by filling 
the rightmost n{0) components with all possible nonnegative fermionic spin 
values, and the rest (if any) with the unique bosonic spin s^^ = 0.) In summary. 
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ii m,n ^ 1 the minimuin energy is given by 



0, m>l, orm = e = l, 

^min ={ I __ . (89) 

max ( A^ — 77,(0), Oj , m = — e = 1 . 



As to the maximum energy, the duahty relation (1731) and the previous equation 
imply that 



n > 1, or n = — e' = 1 , 

^max=<! ,^/ . ^ _ N , . (90) 




A^2_niax(A^-m(0),0) 



n = e 



From Eqs. fl77|) . fl75]) and the last two equations, it easily follows that both 
(^max — lAI'^ ^'^d (/i — £^min)/<7 a^G 0(A^2) as A^ ^ oo, SO that condition (iii) 
above is indeed satisfied. On the other hand, by Eqs. (1771) . (15^ and (MI]) we 
have 



SO that condition (iv) is not satisfied unless 7n = 77, as claimed above. 

We have verified that the approximation (!80|) - (!8T]) is in excellent agreement 
with the numerical data obtained using our exact formulas for the partition 
function for many different values oi m,n ^ 1 and A^ ^ 15. For instance, 
in Fig. [3] we plot the cumulative spacings distribution P{s) versus its ana- 
lytic approximation flHOj) in the cases m = n = —e = e' = l and m = 2, 
n = e = e' = 1 ioT N = 15 and A^ = 30 spins. The respective root mean 
square errors (normalized to the mean) drop from 9.4 x 10~^ and 4.6 x 10"^ 
for A^ = 15, to 2.4 x 10"^ and 1.1 x 10"^ for A^ = 30. It is worth noting that 
the approximation fl5Ul) - flHTl) for the cumulative spacings distribution contains 
no adjustable parameters, since the maximum spacing Smax is completely de- 
termined by Eqs. (1781) . (!89|) and (!90|) . In fact, from the previous equations one 
easily obtains the following asymptotic formula for the maximum spacing in 
the genuinely supersymmetric case m,n'^l: 



3 / ?)Amn — rn? — n^\ 2^1 ^/,. 1 
V27r V (m + 77)"* / 



2 



Thus the behavior of Smax for large A^ is qualitatively the same as in the 
non-supersymmetric cases studied in Ref. [53] . 
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Fig. 3. Cumulative spacings distribution P{s) (blue dots) and its analytic approxi- 



mation (|80p - (j8ip (continuous red line) in the cases m 



n 



-e = e' = 1 (top) and 



m = 2,n = e = e' = l (bottom) for N = 15 (left) and N = 30 (right). 
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A Some useful g-number identities 



In this appendix we shall collect several identities involving g-numbers that 
are used in the simplification of the partition function Z^ of the chain ([5]) 
when either m or n are odd. Given a real number q G (0, 1) and a nonnegative 
integer k, we define the symbols {q)k and [k]g by 



k 



iQ)k = m^-i' 



[k\ 



q^ 



1-? 
with (g)o = 1. The q- factorial of k is then defined as 

k 



k-1 
i=0 



(A.l) 



[kW■ = E[^ 



IQ ■ 



1=1 
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Finally, if k,l are nonnegative integers with k ^ I, we define the q-hinomial 
coefficient 



k 



I 



[kV 



[l],\[k-l],\- 
From the above definitions it immediately follows the useful equality 



(?)a 



{q)M) 



(A.2) 



k-l 



An important identity satisfied by the g-binomial coefficients is the following 
instance of Newton's g-binomial formula [39l Eq. (8)]: 



E 

i=0 



k-l 
q j=0 



(A.3) 



From the previous equation with q replaced by q^ and x by g we obtain the 
formula 

^1 / = n(l + g2*+i). (A.4) 

Replacing q by q^ and setting a; = g~Mn Eq. flA.Sp we obtain the analogous 
relation 



k 

E 

i=0 



E 

i=0 



fc-1 



/-2) = n(i+g'*"': 



(A.5) 



i=0 



The last three identities are used in Section [5] to compute the partition func- 
tionZi:i^\ 



Likewise, if 5q is the g-dilation operator, whose action on a smooth function 
/ : R ^ M is given by {5J){x) = f{qx), we have [39, Eq. (11)] 



k 

E 

«=0 L 



X' = {X + 5qf ■ 1 



-I Q 



where 1 denotes the constant function identically equal to 1. Equating the 
coefficient of x* in both sides of the previous equation we obtain the remarkable 
identity 

^1 =g-(^)r,(l,g,...,g'=-i), (A.6) 

Aq 

where rj(xi, . . . ,Xfc) denotes the elementary symmetric polynomial of degree 
i 'vn k variables, defined by 



Ti{xi,...,Xk) = Yl 



•^ii ■ ■ ■ -^ji 



(A.7) 



In particular, from Eq. flA.6P it easily follows that the g-binomial coefficient 



is a polynomial of degree i{k — i) in q. 
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B Computation of the mean and variance of the energy 



In this appendix we shall compute the mean and variance of the energy of the 
chain i^ for arbitrary values of N, m, n, e and e'. In the purely fermionic case, 
these quantities were computed in Ref [33] using the formulas for the traces 
(of products) of the operators Sij and Si given in Ref [31]. In the genuinely 
supersymmetric case, one should evaluate the traces of products involving 
the more general operators 5*^^ and ST'- Indeed, the traces of S^T and 






(m|n) Q{m,\n) 



"U 



^kl 



which are the only ones appearing in the chain of type A, were 
computed in detail in Ref. [27]. The calculation of the traces involving the 
spin reversal operator S*" is totally analogous, so that we shall limit ourselves 
to listing the final results in Table [BTTI In this table Pm,Pn £ {0, 1} are the 
parities of m and n, respectively, and we have omitted the traces of products 
which reduce to the identity matrix or to one of the entries. 

Table B.l 

Traces of products of the spin operators 5)"* , S^P and 5" . 



Operator Trace 


Sf {m + nf-\epm + e'pr,) 


^(Hn)^ ^Hn) (^_„)(^ + ^)^-2 


SfSf ^rn + n)''-\epra + e'pnf 


S\T^^^SI^\ si^l-^Sf {m-n){m + nr~\ep^ + e'pn) 


giMn)^., siMn)^., (^ + nf-^{ep^ - e'pn) 


Q{ra\n) (j(m\n) n{rn\n) Q{m\n) n{m\n) n{m\n) , „\2f^ , „\N-4: 


Qim\n) ^{m\n) r,(m|n) 5(r?i|n) Mm\n) n{m\n) f^ i „\N-2 



In the first place, calling 



hi 



(e. + e, 



h^=/3C^ 



and using the formulas in Table IB. II for the traces of the spin operators we 
easily obtain 



/^ 



tr7^;:?i") 

{m + n)^ 

m — n 



(m + nV 



T.ih^J+h)+il- '^"'^''^'' )T.h^ 



■i-T^J 



m + n 



The sums in the previous expression were evaluated in Ref. 
result 



(B.l) 
with the 
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Inserting these values into (IB.ip we immediately arrive at Eq. (1771) for the 
mean energy /x. 

The evaluation of the standard deviation of the energy 



tr 



a 



^^(rn|n)y 



(m + n) 



N 



/i 



(B.2) 



which is considerably more involved, is also performed in two steps. We be- 
gin by expanding Eq. (1B.2P and simplifying the resulting expression using 
Eq. (IB.ip and the trace formulas in Table IB.li After a long but straightfor- 
ward calculation we obtain 



a 



[m — n) 



[m + nY 



2E(^?,+^?;)+E^^ 



«7^i 



+ 



+ 



[m + ny \\m + n 



m — n 



{ep„, + e'pnf ( 7 E ^? - E ^if^ij 



«7^i 



'm + riY 






+ 7 : — —2^ [hij + hij)[hik + hik] 



[m + nY 



(B.3) 



«,i,fc 



where the symbol X)' denotes summation over i ^ j ^ k ^ i. The sums in 
the first two lines of Eq. ( IB. 31) , which were evaluated in Ref . [33] , are given by 



2 Y^iK + hi) + Y.hl = -^ (4iV2 + 6iV - 1) 



4^ 



iV 
36 

16 



(B.4) 
(B.5) 



In order to evaluate the last two sums in (]B.3p . we recall from Appendix A of 
Ref. [33j the following identities involving the zeros yi = CF/2 of the Laguerre 
polynomial L^~ and their differences yij = yi — yj : 



^ 1 _ N{N + (3) 
^yf~ P^iP + 1) 
1 _ N 
_ Vi /? ' 



E 



Y,y^ = N{N + ^-l){2N + P-2) 



(B.6a) 

(B.6b) 
(B.6c) 
(B.6d) 
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= ^(Z/.-/?) 


V- yf 

^ 7/2. 


1 ( 2 


^7,2 


4(/5+l) 



-!/f + 2(2iV + /3)a,-/3(/3 + 4) 



yf N{N-l)(^{2N + 5)(3 + 2N + U 



144(/? + 1) 



(B.6e) 
(B.6f) 
(B.6g) 
(B.6h) 

(B.6i) 



Multiplying Eq. flB.6eP by y^ ^ and i/i, summing over i the resulting identities, 
and using flB.6ap - flB.6dl) . one easily obtains 



-,., -W^- gS-^'('^-l)('^ + .''-l)■ (B'^) 



E^ 



Similarly, summing Eq. flB.6fj) over i and using flB.6cP - flB.6dp we get 



(B.8) 



The sum in the third line of Eq. fjB.SP can now be easily computed. Indeed, 



E(^*j + ~hi3)hi = o E ^' ^^ 



«7^i 



-1 



2 •_^^. ViVij 

/3i:i-fi:^=iA'('v-i) 



(B.9) 



where we have used Eqs. (B.6g) and flB.7p . Turning next to the last sum in 
Eq. fIB.Sp . we have 



Y^ {hij + hij){hik + hik) = Y^' 

i,j,k i,j,k 



(vi + yj)iyi + Vk) 



yij ytk 



Y' ^yf Y' ^y^ Y' ^ 
i,j,k yij yik i,j,k yij yif^ i,j,k yij yik 



(B.IO) 



The last sum in the previous equation clearly vanishes on account of the 
identity 



+ 



+ 



1 



{Zl - Z2){Zi - Z^) {Z2 - Z3){Z2 - Zi) (Zg - Zi) (Zg - Z2) 
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valid for arbitrary Zi ^ Z2 ^ z^ ^ Zi. On the other hand, the first two sums 
in the RHS of Eq. (JB.IOP can be computed using Eqs. (lB.6p - (]B.8|) . Indeed, 



%.^ 



4 



E ^ / v^ Vi yj_ 
^,2 I Z^ ^,2 ^,2 



i,j,k yij yik ijLj Vij \ k,k^i "^ik Vij 

= ^Ei(-l/' + 2(2iV + /3)y,-/?(/? + 4))-4^^ 



9 



N{N -l){N-2) 



i^j yij 



where we have used Eqs. (]B.6f|) - (]B.6hp and (IB.Sp . Likewise 



i,j,k yij y^^ i^j y^j v k,k=ii y^^ y^-j , 



^sr^ Vi- ji _sr^ yj_ 



l^ yi-(3 I I _l y, /3 + 1 1 



ijLj yij ^ ijtj yij ^ ijLj yij ^ i^j yij 



on account of Eqs. ( |B.6gD -( lB.6h[) . The LHS of Eq. ( JB.lOp thus reduces to 



Y, {h,, + h,,){h,k + h,k) = -N{N- l)(iV - 2) . 



(B.ii; 



Substituting Eqs. (IRJD . (|R5|) . (IROJ) and (IRTOJ) into (IRSJ) we finally obtain 
the expression (!78|) for a^ in Section [71 
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